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ABSTRACT 

We discuss the effects of a comptonizing corona on the appearance of the reflection 
components, and in particular of the reflection hump, in the X-rays spectra of accret- 
ing black holes. Indeed, in the framework of a thermal corona model, we expect that 
part (or even all, depending on the coronal covering factor) of the reflection features 
should cross the hot plasma, and thus suffer Compton scattering, before being ob- 
served. We have studied in detail the dependence of these effects on the physical (i.e. 
temperature and optical depth) and geometrical (i.e. inclination angle) parameters of 
the corona , concentrating on the slab geometry . Due to the smoothing and shifting 
towards high energies of the comptonized reflection hump, the main effects on the 
emerging spectra appear above 100 keV. We have also investigated the importance 
of such effects on the interpretation of the results obtained with the standard fitting 
procedures. We found that fitting Comptonization models, taking into account comp- 
tonized reflection, by the usual cut-off power law + uncomptonized reflection model, 
may lead to an underestimation of the reflection normalization and an overestimation 
of the high energy cut-off. We discuss and illustrate the importance of these effects 
by analysing recent observational results as those of the galaxy NGC 4258. We also 
find that the comptonizing corona can produce and/or emphasize correlations be- 
tween the reflection features characteristics (like the iron line equivalent width or the 
covering fraction) and the X-ray spectral index similar to those recently reported in 
the literature. We also underline the importance of these effects when dealing with 
accurate spectral fitting of the X-ray background. 

Key words: radiation mechanisms: thermal - X-rays: general 



1 INTRODUCTION 

The presence of secondary components in the spectra 
of Seyfert galaxies and Galactic Black Hole Candidates 
(GBHC), such as an iron line at ~6.4 keV and a reflec- 
tion hump between 10 and 50 keV, superimposed on the 
primary X-ray continuum, is now well established observa- 
tionally ( [Mandra fc Pounds 1994 ). Discovered in the late 



eighties with the G INGA satellite (Matsuoka et al. 1990 



Pounds ct al. 199C), they are signatures of reprocessing of 
the primary X-ray emission in surrounding cold (T < 10 5 K) 
and optically thick matter accreting onto the central engine. 
Indeed, part (~ 10%) of the primary X-ray radiation may 
be Compton reflected at th e gas surface, produc i ng the ob- 
served secondary features (|Bai fc Ramaty 1978 ; Lightman 
fc White 1988|; | George fc Fabian 199l| ; [Matt, Perota fc Piro 
1991; Magdziarz & Zdziarski 1995). The remaining incident 
flux is reprocessed in the UV/Soft X-ray band and is be- 



lieved to form part of the UV-soft X-ray excess (the so- 
called UV bump) generally observed in this class of objects. 

The reflections features are of great importance in test- 
ing theoretical models of the high energy emission in com- 
pact objects since they give crucial (while undirect) con- 
straints on the geometry and on the nature of the emitting 
regions. 



Compton scattering of the soft photons emitted by the 
thick matter on a population of hot thermal electrons is the 
best model to date for the primary high energy emission in 
these objects. Indeed, the so-called corona models, which as- 
sume radiative equilibrium between the thermal comptoniz- 
ing plasma and the underlyimg cold matter (generally an 
accretion disk), can naturally account for the average prop- 
erties of the X-ray emissio n of black holes accretion flows 
(Haardt fc Maraschi 1991). It has also been shown that a 
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patchy geometry, where the corona is disrupted in local- 
ized blobs, appears to be in better agreement with the spec- 
tral variability a nd Luv / Lx r atios observed fo r example in 



Seyfert galaxies (Haardt 1994; Stern et al. 1995). In the case 
of GBHC, like Cyg X-l, simple slab-corona models even fail 
to fit the high signal to noise X-ray spectra of these objects 
and more complex configurations are required (Gierlinski et 
al. 199'||). The real geometry of the corona-disk system is 
thus relatively difficult to constrain in the framework of this 
class of models and the reflection features can give some im- 
portant clues in this respect. 

An alternative model for the high energy emission in 
black hole spectra consists of an advection-dominated ac- 
cretion flow (ADAF) near the central engine. ADAFs are 
expected to occur in low mass accretion rate objects, where 
the hot plama density is low enough to prevent a rapid en- 
ergy transfer from protons to electrons. Most of the accretion 
energy is then advected into the central black hole rather 
than being radiated (Narayan et al. 1998, and references 
therein). 

The amount of reflection is expected to be different in 
these different scenarios. For instance, in a disk-corona con- 
figuration we expect a reflection fraction of the order of 
unity, since the reflected disk, as seen by the hot corona, 
subtends a solid angle of 2n. On the other hand, in the case 
of a truncated disk-ADAF model, the solid angle would be 
smaller due to the lack of reflecting matter in the central 
regions. For the same reason, we expect narrower iron Ka 
fluorescent lines since relativistic effects, which are stronger 
in the vicinity of the black hole, are negligible in the ADAF 
interpretation. 

We thus expect the reflected spectrum to provide a 
strong obervational test to discriminate between these dis- 
tinct models and between different possible geometries. 
However, various effects may complicate the situation. For 
instance, the complex ionization pattern of the reflector can 
strongly mask and/or modify the reflection features in com- 
parison to what we expect in the case of a simple neutral ma- 
terial oreyen single-zone ionization model. Done fc Nayak- 
shin (2001) and Ballantyne, Ross & Fabian (2000) have re- 



cently shown that using single-zone ionization models can 
severely underestimate the reflection normalization if the ac- 
cretion disk is highly ionized. It can even produce an appar- 
ent correlation between the reflection normalization and the 
spectra l pho ton index similar to that reported by Zdziarski 
et al. (1999) in a sample of Seyfert galaxies and galactic 
black hole candidates. The presence of a dynam i c hot corona 
( iReynolds fc Fabian 1997^ |Beloborodov 1999| ; |Malzac, Be- 
loboroc jov fc Poutanen 200( ) will also strongly modify the 



reflection spectrum due to the anisotropic illumination pro- 
duced by relativistic motions inside the corona. 

We present here another important effect which will also 
modify the appearance of the reflection components. Indeed, 
in a corona-disk configuration, we expect that part (or even 
all, depending on the coronal covering fraction) of the re- 
flection features cross the comptonizing plasma before being 
observed. In this case, the reflected photons are also Comp- 
ton scattered in the corona and the shape of the secondary 



spectral components may be significantly modified as com- 
pared to the one expected when no Comptonization in the 
corona is taken into account. 



The effects of the hot corona in modifying the reflec- 
t ion component hav e been already noted by some authors 
(Haardt et al. 1993b). The aim of this paper is thus to study 



more precisely the dependence of these effects on the phys- 
ical (temperature and optical depth) and geometrical (in- 
clination angle) parameters of the hot corona. We are also 
interested in the consequence for spectral analysis of real 
data when these effects are taken into account. We mainly 
focus on the slab geometry where the effects are expected 
to be large since the corona completely covers the reflecting 
material. 

The paper is organized as follows. In section ^| we briefly 
summarize the main characteristics of the Comptonization 
model we used. We discuss the effects of a comptonizing 
corona on the reflection hump shape, varying the corona op- 
tical depth and/or temperature and for different inclination 
angles in section ^. In section ^, we focus on the iron line 
and on how the Comptonization process affects the measure- 
ment of its equivalent width. We also investigate, in section 
H the importance of these effects on the interpretation of 
the results obtained with the standard fitting procedures. 
We discuss and illustrate the importance of these effects in 
section ^| by analysing recent observational results as those 
of the galaxy NGC 4258. We also show that the comptoniz- 
ing corona can produce apparent correlations between the 
reflection features characteristics (like the iron line equiva- 
lent width or the covering fraction) and the X-ray spectral 
index similar to those recently reported in the literature. We 
finaly underline the importance of these effects when deal- 
ing with accurate spectral fitting of the X-ray background 
before concluding. 



2 THE MODEL 

We use the thermal Comptonization code of Haardt (1994, 
hereafter H94). This code computes the angle-dependent 
spectrum of the disk-corona system in a slab configuration 
using an iterative scattering method, where the scattering 
anisotropy is taken into account only up to the first scatter- 
ing order. It also computes, separatly, the reflection compo- 
nents obtained including or not the comptonization effects. 
They will be denoted Rdir and i? C omp, respectively, in the 
following. The continuum will be simply called C. 



The reflection componen t _Rdi r is computed following 
White , Lightman & Zdziarski (1988) and Lightman & White 
( 1988[ ); thus, the spectral shape of the reflected photons is 



averaged over angles. It is well-known that the actual shape 
of the reflection component does depend on the viewing an- 
gle, especially at high energy. Such angular distortion is rel- 
atively complex and is a function of the photon energy. The 
main effects occurs in hard X-rays and soft 7-rays bands, 
where the reflected spectra strongly hardens with increas- 
ing viewing an gle ( Hua fc Lingenfelter 1992 ; Vlagdziarz & 
Zdziarski 1995|). We have checked a posteriori (cf. section 
that the effects produced by such angular 



3.2 



and Fig. 



dependence are always smaller than the one we discuss here. 
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Figure 1. (a) The comptonized (solid lines) and uncomptonized (dashed lines) reflection humps, R CO mp and Rdir, for different values of 
the coronal optical depth, the temperature being fixed to fcT e =50 keV. From bottom to top, t=0.1, 0.3, 0.5, 0.7 and 0.9. (b) Deviations 
(in %) between the outgoing spectra. The solid line and dot-dashed lines correspond to kT e = 50 keV and 350 keV respectively. 



Consequently, in the rest of the paper, we will normalize the 
spectrum at different viewing angles by just multiplying the 
reflection component by the angular function 



following. 



(3m 2 



(3 - 2// + 3^ 4 ) In ( 1 + - 
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neglecting any energy dependence (see Ghisellini, Haardt & 
Matt 1994 for details). It is worth noting that the shape of 
the reflection component averaged over angles is very similar 
to the real shape expected with an inclination angle of 60°. 

The effect of Comptonization on the reflection com- 
ponent will depend mainly on three physical parameters: 
the temperature kT B and the vertical optical depth r of the 
corona and the inclination angle i (or its cosine p = cosi). 
The last two parameters may be combined to give an "ef- 
fective" optical depth t m = r//i, which is the line-of-sight 
optical depth that a photon, emitted at the surface of the 
disk, has to cross in order to escape from the corona with- 
out being comptonized and reach an observer at the viewing 
angle i. 

The soft temperature of the cold matter fcTbb, does not 
play an important role in the problem we study here. In 
fact, in an anisotropic geometry, as the one we are dealing 
with, the spectrum emitted by the corona depends on such 
t emperature because of the presence of an anisotropy break 



( Stern et al. 199?: ; Svensson 1996 ; Haardt, Maraschi fe Ghis^ 
ellini lpT\ |Petrucci et al. 200^ ). Such a break is due to the 
apparent reduced contribution of the first Compton scatter- 
ing order, mainly emitted backward towards the disk, to the 
outgoing flux. The spectrum is then better approximated by 
a convex broken power law where the break energy depends 
mainly on fcT and fcTbb- Changing fcTbb will thus modify 
the shape of the X-ray spectrum impinging on the cold mat- 
ter and consequently the shape of the reflection components 



R diT and R c . 



However, we have checked that these effects 



are negligible in comparison to the ones we are looking at. 
We will thus fix fcTbb to a constant value of 10 eV in the 



3 COMPTONIZATION EFFECTS ON THE 
REFLECTION HUMP SHAPE 

3.1 Dependence on the physical parameters of the 
corona 

In this section we will firstly discuss the effect of the Comp- 
tonization in the corona on the shape of the reflection com- 
ponent varying the optical depth and/or the temperature 
of the corona, but keeping fixed the viewing angle (here 
H — 0.9). According to thermal Comptonization theory, the 
temperature of the corona is directly related to the average 
fractional energy change that a photon undergoes per scat- 
tering, whereas r gives an estimate of the mean number of 
scattering events. Both parameters play thus a major role 
in the Comptonization process. 



3.1.1 Varying r at fixed coronal temperature 

We will first suppose here that the temperature of the corona 
is fixed, equal to 50 keV, while r varies. We have plotted in 
Fig. [j^, the reflection shapes of i?dir and R CO m P (in dashed 
and solid line respectively) for different values of r. 

The increase of the optical depth hardens the X-ray 
primary spectrum, thus modifying the intrinsic shape of the 
reflection component. It thus explains the hardening of the 
uncomptonized reflection 7?<jir between r=0.1 and r=0.9. 
Increasing r also magnifies the effect of the Comptoniza- 
tion. The larger r, the larger the probability of a photon 
to be comptonized. In this process the reflected photons are 
shifted towards higher energies and an increasing deviation 
between iidi r and 7? C om P is seen, at energies below and above 
~ fcT e =50 keV, for increasing r. 

The optical depth of the corona also controls the rel- 
ative intensity of the different Compton scattering orders. 
The larger r, the higher the intensity of the different orders 
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Figure 2. (a) As in Fig. gja, but for a constant Compton parameter y equal to 0.6, i.e. consistent with a slab corona in radiative 
equilibrium above a passive accretion disk. From bottom to top t=0.1, 0.2, 0.35, 0.5 and 0.64 and fcTe=240, 150, 90, 65 and 50 keV. (b) 
As in Fig. [l]b but for a constant Compton parameter y equal to 0.6 (solid line) and 2 (dot-dashed lines). The latter case is consistent 
with a hemispherical corona in radiative equilibrium above a passive accretion disk. 



producing a hardening of the comptonized shape at high en- 
ergy (above ~ 10-30 keV) as observed in Fig. [l]a. 

To estimate quantitatively these effects on the total out- 
going spectrum, we have plotted in Fig. [j]b (solid lines) the 
deviations (in %) between the total spectrum (C + -Rcomp), 
expected if all the reflected photons cross the comptonizing 
corona before being observed, and the total one (C + -Rdir 
predicted when the Comptonization of the reflection hump 
is not taken into account. The most important effects (vari- 
ations of > 50%) occur at very high energies (above 100 
keV). They are due to the hardening of i? C om P for large op- 
tical depths. For the larger optical depth case considered 
(r = 0.9), a factor of ~2 is expected near 1 MeV. At these 
energies, .Rcomp tends to have the same shape as the pri- 
mary continuum and the fractional deviations between the 
two spectra attain a maximum. 



comptonized flux as soft photons which are the seeds for 
Comptonization. In this case, the system can only be in 
equilibrium if th e temperature and optical depth satisfy a 
precise relation ( Haardt fc Maraschi 1991 ). That in turn 
corresponds to roughly constant Compton y parameters, de- 
fined as: 



> »= 4 (S)[ i+4 (S)k^> 



(2) 



which is about 0.6 in a slab geo metry when all the accretion 
p ower is released in the corona ( Haardt & Maraschi 1993a; 



Stern et al. 1995 Svensson 1996) 



In the case of a constant Compton parameter, changes 
in r will be necessarily accompanied by changes in kT e , 
namely a larger optical depth will require a smaller coronal 
temperature. This results in three effects: first, the photon 



index of the primary contin uum correlates with r (Haardt 



The Cdthf)t6tllzati6tl f>r6duces also s'ulallet (< 10%) dif- 
ferences between Rdh and .Rcomp near 10 keV. The larger one 
is still produced in the high optical depth case, as expected. 
It is worth noting that the X-ray detectors are generally well 
sensitive in the 10-20 keV energy range, which can help con- 
straining the degree of Comptonization. 

Also plotted in Fig. [j]b are the results obtained for a 
coronal temperature of 350 keV (dot-dashed lines). We can 



Maraschi & Ghisellini 1997), whereas these two parameters 



arc anti-c orrelated if the temperature is kept fixed (cf. sec- 
tion 3.1.1). Second, the variation of the coronal temperature 



still observe differences between Hair and R c , 



below and 



above ~ kT e . However, such deviations are larger than in the 
previous case, at least at low energy (below 100 keV), since 
the average fractional energy change of a scattering photon 
increases with kT a . These deviations reach ~40% near 20-30 
keV for t = 0.9. Above 100 keV, they roughly behave as in 
the small temperature case. 



3.1.2 Varying t and T c at fixed Compton parameter 

In the case of a disk-corona system, such as the one we 
are dealing with, we expect the comptonizing region and 
the source of soft photons to be coupled, as the optically 
thick disk necessarily reprocesses and re-emits part of the 



modifies the high energy cut-off of the primary continuum, 
which scales roughly as 2kT e . Consequently, it slightly mod- 
ifies the high energy part of the reflection component -Rdir- 
Finally, the change of kT a modifies the mean photon energy 
gain per scattering. 

These different effects produce a hardening of -Rcomp 
at high energy when r decreases as can be seen in Fig. 

This behavior is thus the opposite of what we obtained 
when the temperature is fixed. However, as shown in Fig. 

(solid lines), the deviations between the outgoing spectra 
closely resemble those shown in Fig. [ijb, becoming larger for 
larger optical depth (i.e. smaller temperature) and reaching 
a plateau at high energies (above 200 keV). For comparison, 
we have also plotted the deviation curves obtained in the 
case of a hemispherical geometry in radiative equilibrium 
(dot-dashed lines). In this case, the Compton parameter is 
of the order of 2 (if all the accretion power is released in the 
corona, see Stern et al. 1995). We see that the deviations are 
smaller (by about a factor of 2) than in the slab case. This 
is simply due to the fact that in hemispherical geometry, 
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Figure 3. Plots of R CO mp (solid line), Rdir (dashed line) and R peX rav (dot-dashed line) at two different viewing angles: (a) fi=0.9 and 
(b) ^=0.3. 




Figure 4. (a) The comptonized (solid lines) and uncomptonized (dashed lines) reflection humps, R CO mp and Rdir, f° r different viewing 
angle. From bottom to top /i=0.1, 0.3, 0.5, 0.7 and 0.9. The optical depth and temperature of the corona are fixed to 0.35 and 90 keV 
respectively, (b) Deviations (in %) between the outgoing spectra. In dot-dashes line the case of an hemispherical corona with r=0.35 
and fcT e =200 keV. 



the corona does not cover all the cold disk. Consequently, 
only a part of the reflected photons have to cross the comp- 
tonizing plasma before being observed, reducing the effects 
of Comptonization. 



3.2 Dependence on the inclination angle 

In this section we fix the temperature and optical depth of 
the corona to 90 keV and 0.35 respectively. Those values 
correspond to a Compton parameter y close to 0.6, i.e. con- 
sistent with a slab corona in radiative equilibrium above a 
passive accretion disk (cf. previous section). We focus here 
on the dependence of the Comptonization effects on the sys- 
tem viewing angle. 

As already mentioned in the introduction, the actual 
shape of the reflection hump depends on /i in a relatively 
complex manner which is a function of the photon energy. 
We recall that this angular dependence is neglected in our 
computation of i?<jir, whose shape is constant and corre- 
sponds to the angular averaged one. We have plotted in 
Figs, ^a and ^|b the reflection component Rdir, computed 



by our code (dashed lin e), an d the one computed following 
Magdziarz & Zdziarski ( jlQQgj ) , obtained with the PEXRAV 
model of XSPEC, and noted i? pC xrav (dot-dashed line), for 
two different inclination angles fi — 0.9 (Fig. Ha) and /i = 0.3 



(Fig. |b) 



is computed using angle dependent Green 



functions and takes properly into account the angular distor- 
tion of the Compton reflection. The continuum assumed to 
compute -Rpexrav corresponds to the best cut-off power law 
fit of the spectrum emitted by our slab corona model, that is 
a photon index P ~ 2.2 and a high energy cut-off E c ~ 250 
keV. The primary spectra impinging on the disk are thus 
roughly similar in the two models. We see in Fig. ^a and |^b 
that the main differences between J?dir and Tipoxrav appears 
above 10 keV. -Rpexrav also hardens for increasing i (i.e. de- 
creasing fi). For comparison, we have also plotted on these 
figures the comptonized reflection fi comp computed with our 
model. It is evident that the effects of the Comptonization 
on the reflection hump are generally more important than 
the intrinsic angular ones, especially at high energies. There- 
fore, in the following we will simply use Eq. (|l]) to rescale 
the amplitude of the reflection component as a function of 
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Keeping in mind the limitations of the above approxi- 
mation, we have plotted in Fig. gja, the different shapes of 
Rcomp and 7?dir for different values of the viewing angles. 
Because an increase of i will correspond to an increase of 
the effective optical depth 7>, we expect to observe, for in- 
creasing i, the same effects we observed for increasing r (cf. 
section |3.l| ). However, the main difference is that now the 
shape of the primary spectrum impinging on the accretion 
disk is constant. 

The differences between R C omp and i?dir are relatively 
strong (cf. Fig. ^a). In our case, r M can easily reach values 
of the order of 2 or 3 for large inclination. It thus results 
in large mean photon energy shifts and hardening of R comp 
at high energy. The deviations between Rdir and -Rcomp are 
reported in Fig. ^p. They may be as large as ~30% in the 
10-50 keV range for /i=0.1 whereas they are of the order of 
a few for nearly face-on configuration. 




Optical depth t 

Figure 5. Reduced equivalent width i?Wcomp(/i)/-EWdir(M = 1) 
of the comptonized iron line in function of the optical depth of 
the corona for different inclination angles. From top to bottom fi 
= 1, 0.5 and 0.1 



Above 100 keV the different shapes overlap near 2kT e 
the larger energy that a photon may reach by Comptoniza- 
tion in the corona. It is worth noting that this overlapping is 
a direct consequence of our treatment of the angular depen- 
dence of the reflection hump. Actually, the reflection shapes 
for different inclination angles would be different especially 
at high energies. The deviations would however not strongly 
depend on the real reflection shape and we see that they 
are of the order of 60% near 300 keV for a slab geometry. 
They reach only 30% for an hemispherical corona since the 
Compt onizat ion effects are less important in this case (cf. 
section 3.1.2). 



et al. (1997). We try to quantify a bit more this effect here. 



We note however that the Comptonization effects be- 
yond ~100 keV may not be easily detected. Firstly because 
the high energy instruments are generally not very sensitive 
beyond ~100 keV and thus prevent any good constrains on 
the spectral shape. Secondly, any deviations from the ideal 
case studied here like, for instance, a stratified temperature 
corona, may distort the derived reflected shape, especially 
at high energy, and complicate the spectral analysis. The ef- 
fects at lower energy (near ~20 keV where the instrumental 
sensitivity is very good) would certainly be more relevant 
for any estimation of the comptonization effects. 



4 COMPTONIZATION EFFECTS ON THE 
IRON LINE 

The equivalent width (EW) of the fluorescent Fe Kq line 
produced by cold matter surrounding the hot corona de- 
pends on different parameters, as the elemental abundance 
of the reflecting matter, its inclination angle as seen by the 
observer, and the ionization state of the reflecting surface 
layers (see Fabian et al. 2000 for a review). It also depends 
on the geometry of the corona + cold matter configuration 
(i.e. on the solid angle subtended by the cold matter as seen 
by the X-ray source). The Comptonization of the iron line 
is another process which may aff ect the measurement of this 
EW as pointed by Haardt et al. (19931 ) and discuss by Matt 



Assuming a slab corona geometry above an accretion 
disk of neutral matter and solar abundance, the dependence 
of the EW on /i can be approximated by the simple formula: 



EW a 



»00 



In 2 

= EW dir (iJ,)e- T > 



■/' 



In 1 



(3) 



The term before the exponentia l is a n approximated for- 
mula, given by Ghisellini et al. (1994), of the angular de- 
pendence of the equivalent width due to limb darkening, 
EWdii(l-i = 1) corresponding to the face-on EW of the un- 
comptonized line. The exponential term e~ T>1 gives the prob- 
ability of a line photon to cross the corona and reach an 
observer at the viewing angle i without being comptonized. 
The comptonized ones will be scattered around in the un- 
derlying continuum and will not contribute to the line flux 
anymore. 

We have plotted in Fig. |B| the reduced equivalent width 
EWcomp(fi') I EWdirifi = 1) as a function of the optical depth 
of the corona and for different viewing angles. We see that 
the effects of Comptonization are very strong, especially at 
large inclination angles. For fi — 0.9 the EW is reduced by 
a factor ~3 when r varies between and 1 whereas it be- 
comes rapidly negligible for r > 0.5 if fi — 0.1. We thus 
expect these modifications of the Iron line EW to be a lot 
more observable than the spectral changes of the reflection 
hump detailed in section ^| 



5 APPLICATION TO FITTING PROCEDURES 

In section |§] we have discussed, in an admittedly qualitative 
way, the changes in the observable reflection hump due to 
the presence of a comptonizing corona above the cold re- 
flector (the accretion disc). Here we will try to give a more 
quantitative estimate of such an effect. To do so, we have to 
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rely on some simple and straightforward measure of the re- 
flection component, and compare such a measure in the two 
cases of a bare disc and of a disc engulfed in a hot optically 
thin corona. 

The most widely used model to fit the hard X-ray con- 
tinuum of Seyfert galaxies and GBHC comprises a cut-off 
power law and an uncomptonized reflection (the so-called 
PEXRAV model of XSPEC, Magdziarz & Zdziarski 1995). 
This model is known to be a poor approximation to the true 
Compton continuum since it does not reproduce the (possi- 
ble) anisotropy break (cf. section 2) and may give (mainly 
for high, i.e. larger than 1, optical depth) a poor modelisa- 
tion of the high energy cut-off shape. However, it can help 
to quantify the effects of a comptonizing corona on the ap- 
pearance of the hard X-ray continuum. For this purpose we 
have proceeded in the following way. We have first created 
two sets of simulated spectra using the thermal Comptoniza- 
tion code of H94 in a slab geometry. In the first set we have 
properly taken into account the Comptonization of the sec- 
ondary components in the corona, while in the second one 
we have neglected such effect. Then we have fitted our sim- 
ulated spectra between 1 and 500 keV with the PEXRAV 
model. We have assumed different viewing angles and differ- 
ent values for coronal temperature and optical depth, chosen 
as to keep y — 0.6. For every simulated spectrum, our fit- 
ting procedure gives the best fit value of the spectral index, 
the high energy cut-off and of the reflection fraction. We 
will denote Ri and -B c ,i, respectively, the value of the reflec- 
tion fraction and of the high energy cut-off obtained for the 
comptonized model, and R2 and E Ct 2 the ones obtained for 
the uncomptonized case. The ratio R1/R2 and Ec,2/E c ,i 
can then help to quantitatively estimate the modifications 
of the best fit parameter values due to the comptonization 
effects. 

In Fig. |^ we show the ratio Ri /R2 as a function of the 
optical depth and for different viewing angles i. As expected, 
due to the smoothing effect of the comptonizing corona on 
the reflection shape, Ri is always smaller than R2, because 
the fitting procedure needs smaller reflection normalization 
to fit the comptonized component. Furthermore, their ra- 
tio decreases with increasing i and/or optical depth in the 
corona. This is in agreement with the results of section ^. 

The shifting of the reflection hump -R CO mp towards high 
energies (as discussed in section 3) as a consequence of the 
Comptonization, modifies also the high energy part of the 
spectra and the estimate of the high energy cut-off so that 
it is always larger for the comptonized case. The amount 
of such shift depends on the temperature and optical depth 
of the corona (cf. Fig. |^a) and the higher the temperature, 
the higher the energy of the scattered photons. At the same 
time, however, in our simulations, higher coronal tempera- 
ture corresponding to lower optical depth, the effect of the 
anisotropy break has to be taken into account. When trying 
to fit the curved spectrum with a cut-off power-law, we ob- 
tain a lower value spectral index and, consequently, a lower 
cut-off energy. The net result is shown in Fig. where we 
have plotted the ratio E c ,2/E c ,i of the measured cut-off en- 
ergies as functions of the optical depth for different viewing 
angles. The ratio has a minimum for r ~0.4. At the highest 
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Figure 6. Ratio R1/R2 versus the coronal optical depth for dif- 
ferent viewing angles but keeping the Compton parameter y equal 
to 0.6. From top to bottom /x=0.9, 0.5, 0.2 and 0.1 
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Figure 7. Ratio E c ,2/E c ,i versus the coronal optical depth, for 
different values of fi. From top to bottom fi=0.9, 0.5, 0.2 and 0.1. 

inclination angle (^=0.1-0.2) we are unable to constrain the 
high energy cut-off in the comptonized case E c ,i, below 500 
keV, so we report for the ratio E c ,2/E c ,i=0. 

The main qualitative conclusion of this section is that 
fitting Comptonization models that take properly into ac- 
count the comptonized reflection by the usual cut-off power 
law + uncomptonized reflection leads to an underestimation 
of the reflection normalization and also to an overestimation 
of the high energy cut-off. Quantitatively, these effects will 
of course strongly depend on the quality of the data, espe- 
cially at high energy (above 100 keV). 

6 OBSERVATIONAL CONSEQUENCES 
6.1 The case of NGC 4258 

^From the X-ray point of view, the nearby bright galaxy 
NGC 4258 possesses the general characteristic of the Seyfert 
class, that is: a power law primary continuum with a pho- 
ton index F ~ 2, a relatively large amplitude variability in 
the 3-10 keV band on time scales of a few tens of thousands 
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of seconds and smaller ones (~20%) on time scales of the 
order of an hour (Reynolds et al. 2000, hereafter R00; Fiore 
et al. 2001, hereafter F01) and the presence of a (narrow) 
iron line near 6.4 keV with an equivalent width measured 
by ASCA of 107±37 eV (R00). The BeppoSAX observation 
gives a more poorly constrained value for the EW of 85±65 
eV. Furthermore, following R00, a broad iron line with an 



EW < 200-300 eV may still be consistent with the ASCA 



data. On the other hand, the bad signal to noise above 10-20 
keV has prevented any good detection of a possible reflec- 
tion component. 

NGC 4258 is well known for the detection, from high 
resolution observations, of a water maser. It is expected to be 
produce d by the accreting gas spiraling down to the central 
engi ne ([ciausscn. Hciligman fc Lo 1984; Watson & Wallin 



1994; Mlyoshi et al. 1995). The combination of the central 
mass estimate (deduced from the water maser properties, 
Heernstein et al. 1999; Miyoshi et al. 1995), and the NIR- 
to- X-ray luminosity suggest that NGC 4258 is an AGN in a 
low state with an Eddington accretion rate of 0.0002 (F01). 
The presence of an ADAF seems however to be ruled out by 
the recent BeppoSAX observations, on the base of both the 
measured X-ray spectral shape and X-ray variability, and 
the most natural explanation of the X-ray luminosity of this 
source may be in terms of Comptonization of soft photons 
in a hot corona (F01). The inclination of the disk is large 
and very well constrained i=82±l degrees (Herrnstein et al. 
1999). We thus expect the Comptonization effects on the 
reflection components coming from the disk to be relatively 
important, especially in the case of a slab corona. 

We can thus wonder, from the results obtained in the 
previous sections, what constraints we can put on the ori- 
gin of the reflection features detected in this object. To this 
purpose, we have first try to estimate the physical parame- 
ters (mainly the temperature and optical depth) of the the 
corona. We have thus fit the BeppoSAX data (downloaded 
from the BeppoSAX archive) with H94, assuming an incli- 
nation angle of = 82° (i.e. jj, = 0.14). Since the intensity 
of the hard component is strongly reduced by photoelectric 
absorption below 2-3 keV, we only fit the data above 2.5 
keV. Due to the low signal-to-noise ratio of the spectrum, 
especially at high energy, and the lack of data to constrain 
the soft photon temperature feTbb> r and kT c are badly con- 
strained, r being smaller than 1 at the 90% confidence level 
for 2 parameters, with a best fit value of 0.05. 

These large uncertainties still enable us to make some 
interesting comments. First, let us assume, as considered by 
R00, that the narrow iron line observed in NGC 4258 is from 
the accretion disk. Then it has to originate at relatively large 
radii (R e > 100 Schwarzschild radii) in order to produce the 
observed small line width. In this case, and since the inclina- 
tion angle of the disk is large (we assume that the inclination 
of the inner part of the disk is also of 82°, i.e. the disk is 
not strongly warped) we expect part of the iron line photons 
to be comptonized in the corona. Thus, as shown in section 
4, the EW may be significantly reduced in comparison to 
the one we expect without a comptonizing corona. Let us 
suppose that all iron line photons have to cross the corona 
before being observed (this would be the case if the coronal 



dimension D is large i.e. D > R e as underlined by R00). The 
observed EW of 100 eV would then correspond to an actual 
uncomptonized face-on line EW (the EWdi r (M — 1) param- 
eter in Eq. (|)) larger than 300 eV (~500 eV for r=0.05 
and ~40 keV for r = 0.5 !). Recent estimates of the aver- 
age narrow iron line equivalent width observed in Sey fert 1 
galaxies rather suggest a value of the order ~ 100 eV (Matt 



200C; Lubinski & Zdziarski 2000). This would corresponds 
to an EWdir(M = 1) °f ~200 eV for r=0.5 and a viewing 
angle of 30° (as estimated by Nandra et al. 1997 for Seyfert 
1). These (admittedly qualitative) computations rather sug- 
gest that the narrow iron line in NGC 4258 does not suffer 
Comptonization. The most simple explanation, with the as- 
sumption that the X-ray emitting region has a disk-corona 
configuration, is that the narrow iron line preferentially orig- 
inates from matter not associated with the accretion disk. 

In this case, and following R00, we can also set an upper 
limit on the equivalent width of the uncomptonized broad 
iron line consistent with its non detection by ASCA. R00 
have estimated that the corresponding face-on EW of a 
possible broad line has to be smaller than ~ 200-300 eV to 
be consistent with the data. This estimate takes into account 
the limb darkening and relativistic bending effects but not 
the Comptonization of the iron line photons in the corona. 
If such an effect is also taken into account (and still assum- 
ing an inclination of 82°), the limits are less strict and a 
broad line with a face-on uncomptonized EW~0.8-1.2 keV 
(for r ~0.05-0.1) would still be consistent with the ASCA 
data. Then if there were a broad line present at the level 
seen in most other AGN (~ 100-300 eV, Nandra et al. 1997; 
Lubinski & Zdziarski 2000), we do not expect to see it. The 
data thus do not require that the inner part of the accretion 
disc be truncated as suggested in R00. 



6.2 Correlations between reflection features and 
X-ray slope 

The Comptonization effects we have discussed so far may 
influence the apparent correlation we observe between the 
reflection features characteristics (broad Fe line EW, reflec- 
tion normalization) and the other parameters of the high 
energy continuum. 

In a recent paper LZ00 have analyzed the complete 
available ASCA database to re-examine the issue of the 
strength and width of the Fe Ka line in Seyfert galaxies. 
An interesting results they report is an apparent correlation 
between the broad iron line EW and the spectral index in 
this class of objects. The correlation is actually observed in 
the average spectra of Seyferts galaxies grouped according 
to their spectral hardness. As underlined by the authors, 
this correlation may be naturally explained in the frame- 
work of thermal Comptonization model. Indeed, in this case, 
the power law slope of the X-ray/7-ray spectrum emitted 
by the hot comptonizing corona is related to the rate of 
cooling by incident soft photons. Then the spectral slope 
will be correlated with the broad line EW, provided that 
the main source of the cooling photons is emitted by the 
same medium that is responsible for the observed reflection 
features. Furthermore, this interpretation is consistent with 
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Figure 8. The solid, dashed and dotted-dashed lines are the 90% envelopes of the uncomptonized iron line EW computed for the three 
spectral indexes reported by LZ00, i.e. T =1.62, 1.82 and 1.94 respectively. The hatched area correspond to the intersection of the three 
envelopes. For any fixed value of EWdir J> 100 eV, this area delimits the range of acceptable common values of the coronal temperature 
in the three samples. Alternatively, for any fixed value of kT e , the hatched area corresponds to the acceptable range for the broad line 
equivalent widths. The left and right plots show the case of a slab and spherical geometry respectively. 



the T—R correlation observed by Zdziarski et al. (1999) in 



a sample of Seyfert and galactic black hole candidates (see 
also Gilfanov et al. 1999; Matt 2000). 

We show here that the Comptonization effects in coro- 
nae of different intrinsic properties (i.e. different tempera- 
ture and optical depth) between different objects may also 
contribute to produce this kind of correlation. 
In thermal Comptonization theory, r and kT D are related to 
the X-ra y spectra photon i n dex T through the approximate 



197€ 



Ziziarski 1985) 



rela tion (|Osterbrock 1974 ; Pozdniakov, Sobol & Siuniaev 



r 



(4) 

In .4 y ' 

In this equation, P sc is the average scattering probabil- 
ity, its expression depending on the coronal geometry, and 



A = 1 + 4 



kT c 



+ 16 



t kT c y 
\m c c 2 ) 



is the average photon energy 



,,„ c c 2 \m c c 2 , 
amplification per scattering. In the following, we only treat 

the cases of slab and spherical geometry. The ex pressi on of 

-P sc in these cases can be found in Zdziarski et al. (1994) and 

Wardzinski & Zdziarski (2000) respectively. 



For a given value of the photon index, we can then use 
Eq. (^) to compute the optical depth r for different values of 
kT c . Assuming that the broad iron lines observed by LZ00 
have been comptonized in the corona, we can then compute 
the EW of these lines before Comptonization, EWdir, by 
inverting Eq. ^ (we fix the inclination angle to 45 degrees, 
which is the best fit value obtained by LZ00). In Fig. 
and ^b (corresponding to a slab and a spherical geometry, 
respectively) we have plotted the 90% confidence envelopes 
for EWdir and kT c obtained for each of the three average val- 
ues of the spectral indexes reported by LZ00, i.e. T= 1.62, 
1.82 and 1.94. These envelopes have been computed by tak- 
ing into account the 90% errors of the photon index and 
of the observed EW reported by LZ00. The solid, dashed 
and dotted-dashed envelopes in these figures correspond to 
r=1.62, 1.82 and 1.94, respectively. The hatched surfaces 
correspond to the region of the parameter space where the 



EWdii computed for the different photon index classes are 
equal within the errors. 

Interestingly, the values of EWdir are consistent with 
each other (within the errors) for a relatively large range 
of coronal temperatures and optical depths, particularly for 
the spherical geometry where EWdir — 100-150 eV for kT e J> 
150 keV and r < 0.6. In the slab case, EWdir ^ 100-150 eV 
for 100 <: kT c <i 250 keV and 0.6 ^ r £ 0.1. These results 
thus show that the presence in different objects, all having 
broad lines with similar EW (before Comptonization) , of hot 
coronae with different physical characteristics (i.e. different 
t and kT e ), may result in an apparent EW-T correlation 
(as the one observed by LZ00). 

It is worth noting that the main implication of this in- 
terpretation is that the coronae producing the harder spec- 
tra (smaller F) must also have the larger optical depth. Such 
behaviour necessarely goes with a change of the Compton pa- 
rameter (i.e. a change of the geometrical and/or energetical 
properties of the corona) since for constant y a correlation 
between r and F is expected (Haardt et al. 1997). 

We note that, in this picture, we also expect a corre- 
lation between the observed reflection normalization R and 
the photon index F. Indeed, as shown in section ^, for larger 
optical depth the Comptonization effects smooth the reflec- 
tion hump so that fitting with an uncomptonized reflection 
gives smaller R. This goes in the sense of the results of Z99, 
although, as shown in section 5, we do not expect strong 
effects (<C 30% for inclination angles smaller than 60°) on 
the reflection normalization. 

A possible reason for an anticorrelation between r and T 
may be the following. Suppose that the corona+disk config- 
uration is the one proposed by Z99 i.e. a central hot plasma 
+ cold disc model, the inner radius of the disk being able 
to vary. In this interpretation, the increase of the inner ra- 
dius of the disk, which will produce a decrease of R and 
r, may be due to the evaporation of part of the inner re- 
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gions of the disk (due to some disk instabilities, Meyer et 
al. 2000; Menou et al. 2000; Turolla & Dullemond 2000) in 
the corona, thus resulting in a increase of the corona optical 
depth. Interestingly, in this case, the comptonization effects 
would reduce the apparent reflection component and thus 
could explain the small R values observed at small F, em- 
phasizing the strength of the observed correlation. 

In conclusion, the main point of this section is that, due 
to likely different corona characteristics from one object to 
the other, the comptonization effects may have some influ- 
ence to produce and/or emphasize the observed correlations 
between the reflection features and X-ray slope. 



6.3 Consequences for the X-ray background 

After the most recent surveys (Mushotzky et al. 2000; 
Hasinger et al. 2001) there is now a tantalizing evidence 
that the X-ray Background (XRB) is mainly due to the in- 
tegrated emission of single sources, i.e. AGNs lying at cos- 
mological distances. All the spectral models for the XRB 
assume that the intrinsic AGN spectrum (before any ab- 
sorption in a putative large scale molecular torus) is the 
sum of a power-law continuum and an uncomptonized re- 
flection. We expect that the effects described here can be of 
some relevance for any detailed fitting model for the XRB. 
In Fig. ^|we show, as an illustrative case, an angle averaged 
reflection hump calculated for kT c = 90keV and r = 0.35 
(corresponding to y — 0.6) for both the comptonized (solid 
line) and unco mptonized (dashed line) case. Inde e d, in all 
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current models (IMadau, Ghisellini & Fabian 1994 



Wilman 



k Fabian 1999 ; |Gilli, Salvati fc Hasinger 20010 , for an av- 



erage continuum slope of about 2, the peak of the XRB 
(in vF v ) is associated with the peak of the reflection hump. 
Furthermore, there is evidence that the observed peak of the 
X-ray background is located at slightly higher energies than 
those predicted by the standard models that neglect the ef- 
fect of Comptonization in the corona. 

Although many uncertainties come into the exact de- 
termination of the possible consequences of such effect on 
the XRB shape, as the temperature and optical depth dis- 
tribution in the different sources or the total coronal cov- 
ering fraction, the work we have presented here should in- 
dicate that Comptonization of the reflection component in 
the corona have to be taken into account when dealing with 
accurate spectral fitting of the X-ray background. 



7 CONCLUSION 

In this paper, we have shown that the effects of a comp- 
tonizing corona on the appearance of the reflection com- 
ponents can be relatively important, modifying the shape 
of the reflection hump and the iron line equivalent width 
measurements. We have studied in detail the dependence of 
these effects on the physical (i.e. the temperature kT e and 
optical depth r) and geometrical (i.e. the inclination angle) 
parameters of the corona, mainly focusing on the case of 
a slab geometry. This is an extreme case since the corona 
covers all the reflecting material and consequently the comp- 
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Figure 9. An angle-averaged reflection hump calculated for 
kT c = 90keV and r = 0.35 (corresponding to y = 0.6). The 
solid line represents the comptonized reflection, the dashed the 
direct one. The slight shift of the peak energy of the hump is an 
indication of the predicted effect of the Comptonizing corona on 
the location of the peak energy in the X-ray background. 



tonization effects are large. The results of this study can be 
summarized as follows: 

• Due to the smoothing and shifting towards high ener- 
gies of the comptonized reflection hump, the main effects on 
the emerging spectra appear at energies below and above 
~ kT c . 

• They are larger for larger optical depth of the corona 
and/or larger inclination angle of the corona-disk configu- 
ration. 

• Fitting Comptonization models taking into account 
comptonized reflection by the usual cut-off power law + 
uncomptonized reflection models, leads to an underestima- 
tion of the reflection normalization and an overestimation of 
the high energy cut-off. 

• The Comptonization effects may strongly reduce the 
equivalent width of the observed iron line, especially at large 
inclination angles. 

These effects may have important consequences on the 
physical interpretation of the presence and/or absence of re- 
flection features in astrophysical objects. As an example, we 
have studied the case of the galaxy NGC 4258. With the as- 
sumption that the X-ray emitting region has a disk-corona 
configuration and given the high inclination angle of the ac- 
cretion disk in this source, the comptonization effects enable 
to find less strict limits on the EW of a possible broad Iron 
line, then explaining its non-detection. Besides, in this pic- 
ture, the narrow iron line observed in this object is believed 
to originate preferentially from matter not associated with 
the accretion disk since, in the contrary case, it would re- 
quire a primary (i.e. before Comptonization) iron line with 
a (unlikely) too large EW. 

We also find that the presence of a comptonizing corona 
can produce and/or emphasize the correlations between the 
reflection features (like the iron line equivalent width or the 
covering fraction) and the X-ray spectral index. Then, sim- 
ilarly to the effects produced by a bad modelisation of the 
complex ionization pattern expected at the surface of the X- 
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ray-irradiated reflecting material, these different effects are 
to be properly taken into account to correctly interpretate 
these correlations. 



Finally wc underline the pnaaihlo itnpnrtanrc nf tVio 



Comptonization effects on the reflection shape when deal- 
ing with accurate spectral fitting of the X-ray background. 
Indeed, there is evidence that the observed peak of the X-ray 
background is located at slightly higher energies than those 
predicted by the standard models that neglect the effect of 
Comptonization in the corona. 
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